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Abstract 

Pc: i .ron- electron pair radiation is examined as a mechanism that could be 
responsible for the impulsive phase emission of the March 5, 1979 transient. 
Synchrotron cooling and subsequent annihilation of the pairs can account for the 
energy spectrum, the very high brightness, and the ~ 0.4 MeV feature observed from 
this transient, whose source is likely to be a neutron star in the supernova 
remnant N49 in the Large Magellanic Cloud. In this model, the observed radia- 
tion is produced in the skin layer of a hot, radiation-dominated pair atmosphere, 
probably confined to the vicinity of the neutron star by a strong magnetic field. 

The width of this layer is only about O.linm. In this layer, ~ 10^ generations of 
pairs are formed (by photon-photon collisions) , cooled and annihilated during the 
~ 0.15 sec duration of the impulsive phase. The very large burst energy implied 
by the distance of the LMC, and its very rapid release, are unsolved problems. 

We mention, nonetheless, the possibility of neutron star vibrations, which could 
transport the energy coherently to the surface, heat the atmosphere mechanically 
to a hot, pair-producing temperature, and have a characteristic damping time roughly 
equal to the duration of the impulsive phase. 
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I. Introduction 

Gamma ray lines have recently been detected from at least two gamma- ray 
bursts. Line emission identified as redshifted 0.511 MeV positron annihilation 
radiation was seen in the March 5, 1979 event (Mazets et al., 1979; Mazets and 
Golenetskii, 1979) , and two lines, one identified as redshifted 0.511 MeV emission 
and the other as redshifted 0.847 MeV Fe^ deexcitation radiation, were seen in 
the November 19, 1978 burst (Teegarden and Cline, 1980), The redshifted 0.511 MeV 
feature in the November 19, 1978 burst was also seen by Mazets and Golenetskii 
(1979). These observations have important implications for the origins of gamma ray 
transients. In this paper we shall focus on the March 5, 1979 transient, because 
of the special problems and interesting physics associated with the probable iden- 
tification of the source of this outburst with a neutron star in a supernova remnant 
in the Large Magellanic Cloud (Evans et al., 1980). In another paper, Teegarden 
and Cline (1980) discuss some of the implications of their observation of red- 
shifted line emission in the November 19, 1978 burst. 

Cosmic gamma ray line emission has also been observed from a variety of other 
sources: solar flares (Chupp et al. , 1973; Hudson et al., 1980), the galactic 

center region (Leventhal, MacCallum and Stang, 1978), and possibly from an as-yet- 
undefined class of transient phenomena (Jacobson et al. , 1978). These observations, 
and their theoretical interpretations, have recently been reviewed in considerable 
detail (Lingenf elter , Higdon and Ramaty, 1978; Ramaty and Lingenfelter , 1979; 

Ramaty, Kozlovsky and Lingenfelter, 1979; Vedrenne, 1980). 

The March 5 transient was observed by instruments on 9 different spacecraft 
which form an interplanetary burst network (Mazets et al. , 1979; Barat et al. , 1979; 
Cline et al. , 1980; Evans et al. , 1980; see also Cline, 1980 for a review). The 
position of the burst, determined by this network to an accuracy of an arc minute, 
coincides with the supernova remnant N49 in the LMC. This is the first identifi- 
cation of a gamma ray transient with any known astronomical object. 
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The March 5 transient has several unique characteristics that set It apart 

from the general class of gamma ray bursts: 

-3 2 

The peak flux of ^ 10 erg/ cm sec above 0.03 MeV (Mazets and Golenetskii, 

1979), was an order of magnitude greater than that of any other gamma ray burst. 

At the distance of the LMC (55 kpc) this intensity implies a luminosity of 
44 

~ 3x10 erg/sec. 

-4 

The exceedingly short rise time of this impulsive phase, < 2x10 sec (Cline 
et al., 1980), suggests that the size of the emission region does not exceed 
~ 75 km, the light travel distance for this time. This requires that the source 
be a compact object. 

The duration of the impulsive phase, ~ .15 sec (Cline et al., 1980) was also 

much shorter than that of typical gamma ray bursts which last for several seconds. 

~4 2 

The energy fluence in the impulsive phase was ~ 2x10 erg/cm (Mazets and 

Golenetskii, 1979) which at the distance of the LMC gives an impulsive energy 
43 

release of ~ 7x10 erg. 

The impulsive phase was followed for at least three minutes by much less intense 

pulsed emission with a period of 8.0 +0.05 sec (Cline et al., 1980). The average pulsed 
-5 2 

flux of ~ 10 erg/ cm sec decreased exponentially with characteristic time of 

~ 50 sec (Cline et al. , 1980), so that the fluence in the pulsed phase was about 
twice that in the impulsive phase. The existence of these pulses suggest emission 
from the polar caps of a rotating neutron star. 

In addition, the March 5 transient was followed by three other outbursts , apparently 
from the same source direction, on March 6, on April 4, and on April 24. The peak 
Intensities of these outbursts, however, were much lower than that of the March 5 
event, and decreased with each burst (Mazets and Golenetskii, 1979) . 

Finally, the energy spectrum of the March 5 transient was also quite different 
from that of other bursts. Mazets et al. (1979) have measure' 1 the spectrum for the 
first 4 seconds, an interval dominated by the impulsive phase, and for tne subsequent 
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28 seconds of pulsed emission. The measured spectrum of the first 4 seconds 
is shown in Figure 1. The low energy (< 0.3 MeV) part of this spectrum could 
be fitted by an exponential with a characteristic energy of 0.03 MeV. Above 
~ 0.3 MeV, however, the burst spectrum hardens, showing a broad spectral feature 
peaked at about 0.4 MeV. This feature could be gravitationally redshifted 0.511 MeV 
annihilation radiation produced close to the surface of a ~ 1M^ neutron star with 
a surface redshift z 0.28. 

The gravitational redshift, together with the rapid rise, the 8 sec pulsation, 
the pulse structure, and the association with a supernova remnant strongly suggest 
that the source of the March 5 transient was indeed a neutron star. 

Despite the fact that the burst source position coincides with that of the 

supernova remnant N49, there are several apparent problems, if the source is 

located in the LMC. These concern the origin and rapid release of the very large 

amount of energy emitted in the transient, and the very high efficiency of the 

radiation mechanism. The required mechanism appears to violate the black-body 

44 

limit in that it must provide a luminosity of several times 10 ergs /sec at 
relatively low photon energies (tens of keV) from a small emitting area (a neutron 
star surface). Indeed, several authors (Mazets et al. , 1979; Helfand and Long, 1979; 
Aharonian and Ozernoy, 1979; Petschek and Colvin, 1980) have argued that the source 
cannot be in the LMC. Mazets et al. (1979), assuming that the emission results from 
accretion on a neutron star, argue that the Eddington limit on the luminosity requires 
that the source must be closer than ~ 300 pc. Aharonian and Ozernoy (1979) applying 
Shapiro and Salpeter’s (1975) spherical accretion model, concluded that the source 
was most likely only 5 to 20 pc away. Helfand and Long (1979) argue that if the 
continuum radiation is due to thermal bremsstrahlung and the teature at ~ 0.4 MeV is 
due to e + -e~ annihilation, then for the emission region to be optically thin to the 
annihilation radiation the source would have to lie within 1 pc. But even if these 
are not the relevant energy sources or emission mechanisms, Aharonian and Ozernoy 
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(1979) also argue that the observation of MeV gamma rays in the burst limits 

its distance to less than 200 pc, because of the opacity of the source due to e + -e~ 

pair production. This argument is essentially the same as that of Cavallo and Rees 

(1978) and Schmidt (1978) for bursts in general. In addition to difficulties with accretion 

models, Petschek and Colvin (1980) point out that even for sources internal to 

the star, the energy could not diffuse to the star's surface in a short enough time. 

In che present paper we do take the point of view that the source of the March 
5 transient is in the LMC, based on the rather convincing positional identification 
of the burst direction with N49. We find that the arguments against an LMC origin 
based on radiation mechanisms are all removed in our e+-e synchrotron cooling and 
annihilation model. Regarding the energy source, its very large magnitude 
requires that it be interior to the star, most likely of gravitational origin. 

The release of this energy clearly cannot be by photon diffusion. An interesting, 
although as-yet not well researched possibility for the March 5 transient, is neutron 
star vibration which could, in a coherent and rapid fashion, communicate the 
internal energy to the upper crust and atmosphere. Decay of these vibrations by 
gravitational radiation could, in fact, account for the duration of the impulsive 
phase (D. Kazanas and P. Meszaros, private communication 1980). 

In the present paper, however, we shall deal mostly with the radiation mechanism. 

In Section II, we provide a qualitative description of the model, while Section III 
provides a quantitative discussion of the radiating layer; for we find that essentially 
all the observed emission should be produced in a thin layer of unit optical depth 
in which synchrotron cooling of e^-e pairs converts hot (~ MeV) radiation from the 
lower lying optically thick regions into the observed emission. We summarize our 
discussion in Section IV where we also discuss some additional difficulties and 
problems of the model. 


• MiG INAL PAGE IS 
; ;P POOR QUALITY 
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II. Qualitative Description of the Model 

A variety of sources of energy for gamna ray bursts from neutron stars 

have been proposed (see review by Ruderman 1975), but most of these appear to 

have insufficient energy to account for the March 5 transient at the distance of 

the LMC. It is beyond the scope of the present paper to discuss in detail 

possible energy sources for this transient. We merely mention the possibility 

(S. Bonazzola, private communication 1980) of a phase transition in the neutron 

star interior (e.g. Hartle, Sawyer and Scalapino, 1975) which could release an amount 

of gravitational energy that is a significant fraction of the neutron star’s binding 

energy. As mentioned above, the energy could be stored in neutron star vibration 
-4 -3 

(period 10 to 10 sec) whose dissipation by gravitational radiation (Thorne, 1969; 

Detweiler, 1975) could account for the ~ 0.15 sec decay time of the impulsive phase 

(D. Kazanas and P. Meszaros, private communication 1980). In such a scenario, 

mechanical energy can be deposited in the upper layers of the neutron star crust 

and atmosphere leading to hot radiation capable of producing positron-electron 

pairs. As we shall see below, cooling and annihilation of such pairs should be 

responsible for essentially all the gamma- ray and hard x^ray emission of the 

transient without encountering the difficulties pointed out in the Introduction. 

Regarding the radiation mechanism, the most obvious difficulty appears to be 

that, while on the one hand, to account for the observed bolometric luminosity 
44 

(L 3x10 erg/sec) a temperature kT > 0.06 MeV is required (at this kT a neutron 
star of 15km radius could emit blackbody radiation at the required L) , on the other, 
the low energy part of the observed spectrum (<0.3 MeV, see Figure 1) would onJ y 
be consistent with black body emission if kT were less that about 0.03 MeV. More- 
over this temperature would be too low to account for the high energy par^ of the 
spectrum, and it would also be too low to produce the e + -e pai^s that appear to be 

responsible for the 0.4 MeV feature. The pairs might be produced in a separate 

2 

spatial region of higher temperature (kT ^ me ), but in this case their 
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annihilation emission would be inconsistent with the data, because its width 
would greatly exceed the observed width (< 0.3 MeV). 

In the model that we propose, synchrotron cooling of hot (initially MeV) 
positron-electron pairs and their subsequent annihilation are responsible for 
the entire observed spectrum of the impulsive phase. Hot radiation, possibly 
resulting from mechanical heating due to neutron star vibrations, produces 
e+-e~ pairs. In the absence of cooling, such pairs annihilate to produce new 
photons of the same temperature, which in turn produce new pairs and so on. 

However, in the presence of a strong magnetic field, synchrotron cooling can be 
faster than annihilation. As we shall see below, for this to occur, B has to 
be about ]0^ gauss. Since in such a magnetic field electrons and positrons of 
MeV energies produce radiation of tens of keV, the cooling effectively transforms 
the hot pair-producing radiation into cooler radiation devoid of pairs. This 
transition must take place in a layer whose thickness does not exceed the photon- 
photon pair production mean free path, because the annihilation radiation from the 
pairs that have already cooled cannot produce new pairs. Since the Compton cross 
section at energies from tens of keV to an MeV is comparable to the maximum yy 
pair production cross section, and since (see below) synchrotron self absorption 
begins to affect the spectrum only below ~ 0.03 MeV, the opacity of the transition 
layer is about unity, and hence essentially all the observed radiation should 
originate in this layer. More specifically, the pairs produced in the transition 
layer by hot photons attempting to emerge from lower lving levels produce the low 
energy 0.3 MeV) part of the observed spectrum by synchrotron radiation. Having 
lost the bulk of their initial kinetic energy by synchrotron cooling, the pairs 
establish an equilibrium temperature 0.1 MeV) where the energy loss to the 
magnetic field is balanced by Comptonization. Annihilation at this temperature 
produces the high energy part of the spectrum. 

This process of synchrotron cooling and subsequent annihilation of e + -e pairs 
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overcomes the apparent difficulties of the radiation mechanism as follows: 

Since the energies of the radiating particles (positrons and electrons of 
MeV energies) are high, this mechanism can produce the observed luminosity 
by emitting at a rate well below the black-body limit. The radiation intensity 
in the transition layer is in fact much lower than the black body value for MeV 
temperatures, thereby allowing the nonequilibrium cooling which is an essential 
ingredient of the model. 

The shape of the low energy part of the spectrum is not representative of 
the temperature of the radiation. Rather, it is synchrotron radiation of much 
higher energy positrons and electrons produced in the transition layer by hot 
radiation from lower levels. The photons from these levels do not contribute 
significantly to the observed spectrum because they are absorbed in the transition 
layer. 

The high temperature (~ MeV) just below the transition layer is entirely 
sufficient to produce copious numbers of e + -e pairs. But the Doppler broadening 
of the annihilation feature is greatly reduced because the pairs are cooled before 
they annihilate. This cooling, of course, produces the low energy radiation which 
constitutes the bulk of the observed luminosity. 

The model also overcomes the other difficulties that hive been mentioned by 
previous authors. Obviously, synchrotron radiation in a strong magnetic field can be 
a much more efficient radiation mechanism than bremsstrahlung. Regarding the 
YY opacity, photon-photon pair production would render the source opaque to high 
energy photons (greater than several MeV), since these would interact with the hard 
x rays above the transition layer. But there is not much evidence in the data 
for photons above ~ 1 MeV (Mazets et al. 1979), ind we would, in fact, not expect 
such photons if the source of the March 5 transient is in the LMC. 
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III. Parameters of the Transition Layer 

The transition layer, which in our model produces essentially all of the 

observed radiation, should be strongly radiation dominated, i.e. the density 

of pairs produced by photon-photon collisions in the layer should greatly exceed 

that of the ambient matter there. Because ambient matter at an MeV has a large 

scale height (~ lO^ctn for H around a 1M^ neutron star) its density must be 

20 -3 

relatively low 10 cm ) in order that the region be transparent. Where it 

not for synchrotron cooling and annihilation, an e+-e plasma would have an even 

larger scale height at this temperature. But the density of the pairs should in 

fact drop abruptly at the transition layer because there is almost no further 

pair production above it. As we shall see below, the width of this layer is only 

26 — 3 

a fraction of a um leading to a pair density in the layer of about 2x10 cm . We 
have obtained these estimates as follows: 

The annihilation of e + -e~ pairs in the transition layer produces the high 
energy (> 0.3 Me V) part of the observed impulsive phase. The densities n+ and n_ 
should therefore satisfy 

a n 2 Ah (r/d) 2 - F a , (1) 

where n n+ n_, a is the annihilation rate coefficient, Ah is the thickness 

2 

of the layer, and F a 250 photons/cm sec is the observed photon flux above 0.3 MeV 

in the impulsive phase attributed to annihilation radiation. At relativistic 

temperatures u is temperature dependent. But since the pair temperature just 

before annihilation should not exceed about 0.1 MeV (otherwise the width of the feature 

-15 2 

would exceed that observed), cr is essentially constant, a** 7.5x10 cm sec (e.g. 

Bussard, Ramaty and Drachman 1979). Since Compton scattering should not smear 

out the line, the thickness of the layer should be about (2nac) \ where in the 

-25 2 

energy ran^e of interest a c is about 3x10 cm . For a neutron star radius 
r * 1.5xl0^cm and d » 55 kpc, eq. (1) yields n *• 2.5xl0 2 ^cm ^ and Ah 10 2 cm. 
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The initial kinetic energy of the pairs, y 0 mc , ae they are produced 
in the transition layer, can be estimated from the data of Figure 1. Since 
in our model the low energy part of the observed energy flux density (E <0.3 MeV) 
is due to the loss of this kinetic energy, we should have, approximately, 

me 2 (y -1) - (F.)" 1 f 0,3 M#V dE E F(E) . (2) 

0.03 MeV 

The data of Figure 1 then yields y Q »• 7. An initial particle energy of about 

3 MeV indicates that the radiation Just below the transition layer should have 

a temperature of this magnitude. The density of this radiation, however, is much 

lower than the black body density at MeV energies. This can be seen by comparing 

26 3 

the pair density deduced above (~ 2.5x10 cm ) with the equilibrium pair density 

at such temperatures (e.g. Landau and Llfshitz, 1958) which turns out to be about 

5 orders of magnitude higher. For this reason, rapid cooling in the outermost layer, 

for example by synchrotron radiation, should be possible. 

The magnetic field B in the layer is estimated to be about 10 11 gauss from 

the requirement that electrons and positrons with y »• 7 produce synchrotron emission 

with photon energies characteristic of r v j low energy part of the observed spectrum. 

To cool the pairs from an initial y •** 7 to y < 1.2 (kT <0.1 MeV), the synchrotron 

o 

loss time, t g — dy/(dy/dt) s , must be shorter than the annihilation time, 

.i Y 

* ln the case of strong pitch angle scattering (which wc assume but 
which can also be expected in view of the very strong turbulence that accompanies 
the vibrational heating) , 

(^)g - 2xl0' 9 B 2 (y 2 -1) (sec -1 ). (3) 

-14 

With the above values of y 0 , y and B f eq. (3) yields t # 5x10 sec which is 

26 - 3 

about an order of magnitude shorter than t a for n 2.5x10 cm . Thus, the pairs 
could indeed annihilate after losing the bulk of their kinetic energy. 
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Th« synchrotron photon flux density produced by the cooling electrons 
and positrons in the transition layer can be evaluated from 

F g (E) - F a E" 1 p(E, Y )/(d Y /dt) eff , (4) 

where (dy/dt) e ^ is the effective* energy loss rate of the particles (synchro- 
tron loss minus possible energy gains), and p(E,y) is the synchrotron emissivity 
per particle per unit pheton energy (e.g, Ginzburg and Syrovatskii, 1964). Eq. (4) 
is insensitive to the lower limit of the integration since most the synchrotron 
emission is produced by pairs close to their initial energy. Possible energy 
gain mechanisms that could affect (dy/dt)^^ are synchrotron selfabsorption and 
Coiitptonization. 

We have evaluated the sclfabsorption opacity for the piesc'ot parameters, 
and we find that the transition layer becomes opaque at photon energies below about 
0.025 MeV. From eq. (4) we also find that if (dy/dtJ^j. - (dy/dt)^, about half 
the synchrotron energy loss is in photons below 0.025 MeV, Therefore, a significant 
fraction of this energy is reabsorbed, and this amounts to an effective energy 
gain for the particles This effect can be taken into account in an approximate 
manner by setting (dy/dt) ££ equal to about half (dy/dt) s * where the latter is 
given by eq. (3) . 

Comptotlization of the pairs is caused by the hot (~ 3 MeV) photons which also 
produce the pairs in the layer. The heating rate by these pnotons is approximately 
given by 

(mc2) l * (5) 
44 

where L 3x10 erg/sec is the luminosity of the impulsive phase For the 
present parameters, (dy/dt)^ is much less than the synchrotron loss rate at MeV 
energies, approximately equals 0.5 (dy/dt)^ fer y 1.2, and exceeds it at lower 
energies. The equilibrium temperature, kT, of the pairs in the layer should 
therefore be of the order 0.1 MeV. But since eq. (4) is sensitive mostly to 
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higher pair energies. Comp ton! zat ion has essentially no effect on (dy/dt) 
as applied to eq. (4). 

The solid curve in Figure 1 has been calculated from the expression 

F(E) - (l+z)[F (E / ) + F (n mc 2 kT)“'2~ exp[-(E' -me 2 ) 2 / (kT me 3 )]} x (4/0.15), (6) 

where the observed photon energy, E, and the emitted energy E' are related 

by E' = (1+z) E with z = 0.26. The first term, F g (E), is from eq. (4) with 

11 2 
y o = 7, B = 10 gauss, dy/dt =0.5 (dy/dt) g and ■ 250 photons/cra The 

second term assumes a gaussian form for the annihilation feature, justified by 
the small fraction of annihilations in flight. The width of the feature is 
determined by kT, which, as discussed above, reflects the balance between synchro- 
tron loss and Comptonization. The curve of Figure 1 is based on kT = 0.075 MeV. 

26—3 

(We note that a pair density of 2.5x10 cm is well below the degeneracy limit 
at this temperature). The factor 4 sec/0.15 sec transforms the 4 sec measurements 
of Mazets et al. (1979) into the expected flux density during the 0.15 sec impulsive 
phase, under the assumption that the contribution of the pulsed phase during 
the first 0.15 seconds is small. 

As can be seen from Figure 1, the solid curve does provide a reasonable fit 
to the data, m particular in view of the many uncertainties of the problem and the 
extreme complexity of the physical system under consideration. The excess flux 
density at the high energy end of the spectrum could be due to those hot photons 
that are not stopped in the transition layer, while the excess at the low energies 
could be the contribution of the pulsed phase. 
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IV. Summary and Conclusions 

We haw provided what amounts to a mostly qualitative discussion for the 
origin of the March 5, 1979 transient. We take at face value the identification 
with the supernova remnant N49 in the Large Magellanic Cloud, and we assume that 
the burst is produced by a catastrophic event within a neutron star such as a 
phase transition in its interior. This should lead to the release of a large 
atrount of gravitational energy which, nonetheless, would still be much smaller 
that the binding energy of the star. The transport of this energy to the surface 
is an unsolved problem, although it appears that neutron star vibrations 
are a viable possibility: they could not only provide a readily available, 

coherent form of energy storage, but also heat the atmosphere mechanically to 
the MeV temperatures that are necessary to provide both the high efficiency of 
the radiation mechanism and the e + -e pairs that seem to be responsible for the 
~0.4 MeV spectral feature. Decay of these vibrations by gravitational radiation 
could explain the duration (~ 0.15 sec) of the impulsive phase of the burst. 

The radiation mechanism is extremely efficient in that it produces a very 
high luminosity in relatively low energy photons, and this appears to violate 
the black body limit. We overcome the difficulty by assuming that the 
radiation is nonthermal emission from high energy MeV) pairs, but that these pairs 
are cooled by a strong magnetic field before they annihilate. The observed 
spectrum, therefore, is a combination of synchrotron and annihilation radiations. 

The magnetic field may have, in addition, another very important role: it could 

help to confine the ambient atmosphere of the star fi om being blown away by the 
super-Eddingtonian luminosity of the impulsive phase. 

We have provided no treatment of the pulsed phase or of the subsequent outbursts 
on March 6, April 4 and April 24. The pulsed phase could be due to black body 
radiation at kT ~ 30 keV from the polar caps of the star, and the subsequent events 

could be smaller aftershocks of the initial event. But we defer any further 
discussion of these problems to future papers. 
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Figure Caption 


Energy spectrum of the March 5 transient during the first 4 seconds 
(Mazets et al. 1979). The curve is the calculated spectrum resulting from 
the sychrotron cooling and subsequent annihilation of e + -e pairs injected 
with initial kinetic energy of 3 MeV into a strong but disordered magnetic 
field of 10 11 gauss. The pair temperature after cooling is 0.075 MeV. 
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